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We derive the representation of the two-spinon wavefunction for the Haldane-Shastry model in 
terms of the spinon coordinates. This result allows us to rigorously analyze spinon interaction and 
its physical effects. We show that spinons attract one another. The attraction gets stronger as the 
size of the system is increased and, in the thermodynamic limit, determines the power law with 
which the susceptibility diverges. 
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Interacting spin-1/2 antiferromagnetic spin chains in 
1 dimension exhibit low-lying excitations carrying spin- 
1/2, called spinons §. The Brillouin zone for one spinon 
IS halved and spinons are semions, i.e., particles 

vifith statistics half that of regular fermions [^1). The 
large-scale physics of a generic 1-d antiferromagnet with 
short-range interaction is given by a spinon gas The 
corresponding energy for an A^-spinon solution is the sum 
of the energies of each isolated spinon, plus corrections 
that go to zero in the thermodynamic limit. 

The additivity of the energy is usually claimed as an 
evidence for a spinon gas to be an ensemble of free 
semions In this letter we challenge this idea by 

carefully analyzing the interaction between spinons in an 
exact solution of a particular 1-d antiferromagnet: the 
Haldane-Shastry model (HSM). The HSM is a system of 
spins on a circular lattice interacting via an antiferromag- 
netic interaction inversely proportional to the square of 
the chord between the corresponding sites. The Hamil- 
tonian is given by 



Hhs = J(^) 



Z/3\ 



(1) 



where Za—e'xp{2Tria / N) and a is the lattice site The 
HSM is the simplest exactly-solvable interacting antifer- 
romagnet in 1-d. It is the prototype of a 1-d spinon gas 
since it does not take marginal logarithmic corrections, 
in contrast, for instance, with the behavior of the Heisen- 
berg model 

Many-spinon solutions of the HSM have been con- 
structed 1^ in analogy to the corresponding spinless con- 
tinuum version of the model However, from the 
corresponding "plane wave" representation of the many- 
spinon wavefunction, the persistence of a spinon interac- 
tion in the thermodynamic limit is not at all transparent. 
Indeed, in the thermodynamic limit the energy is the sum 
of the energies of each isolated spinon. However, the in- 
teraction between spinons is "hidden" in the nontrivial 
relation between the canonical momenta, which label the 
states, and the kinetic momenta, which determine the 
energy [|. 



In this paper we work out the real-space coordinate 
representation for two-spinon eigenstates of Tins a-nd the 
corresponding Schrodinger equation. Spinon interaction 
and its nature follow straightforwardly from the behav- 
ior of the exact solution of this equation. In Figj^ we 
plot the result. While at large separations the probabil- 
ity amplitude is independent of spinon separation, as it 
is appropriate for noninteracting particles, at short sep- 
arations there is a huge enhancement. Such a resonant 
enhancement is a clear evidence for a short range, attrac- 
tive interaction between spinons. As we show in Fig||, 
this enhancement gets sharpened as the number of sites 
increases, at odds with the belief that spinon interaction 
and its effects disappear in the thermodynamic limit. 




FIG. 1. Square of two-spinon wavefunction |pmn(2)| de- 
fined by Eq. @ for the case of = 300, m = N/2 - 1, and 
n = 0. At large separations the probability oscillates between 
and 2 and averages to 1. The inset shows this function close 
to the origin for N = 200, 400, and 600. The value at the 
origin diverges in the thermodynamic limit. 

Spinon dynamics determine the low-energy physics of 
the HSM. 1-d interacting antiferromagnets do not order 
and, accordingly, the spin-1 spin-wave (SW) is an unsta- 
ble excitation of the HSM. The SW is absolutely unstable 
at any energy and momentum against decay into a spinon 
pair ||l|. This causes non-analyticities in the SW prop- 
agator, the dynamical spin susceptibility (DSS) Xqi^)- 
Xqi^) develops a branch cut at the threshold energy for a 
SW and a broad continuum above this threshold. Broad 
spectra have been observed by means of neutron scat- 
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tering on quasi 1-d samples |L(]|] , which experimentally 
substantiates this scenario. However, the continuum is 
not flat, as would be the case is it were a spinon joint 
density of states, but rather has a divergent square root 
edge. We shall show that it is the spinon interaction 
which makes the matrix element for the decay of the spin 
wave into spinon pairs huge at threshold, and causes this 
divergence. We explicitly prove that, in the thermody- 
namic limit, the spinon attraction turns into the square 
root divergence in the DSS. Spinon interaction and its 
relation to the DSS are the main result of our work. 

Let us begin with some basic results from the HSM. 
In the even-iV case the ground state of Hhs (eq.(0)) is a 
disordered spin singlet, whose wavefunction is given by 



M 



M 



*GS(Z1, ■ • ■ , Zm) = - ^j)^ 



(2) 



i<j 



where M = N/2 and the {j}'s denote the positions of |- 
spins, all the others being |. The corresponding energy is 
given by Eos = - J(7rV24)(iV + 5/iV) Elemen- 
tary excitations above 'ilcs Siie spinons-spin-1/2 defects 
in the otherwise featureless disordered sea. A | spinon 
localized at a can be thought of as a singlet sea where 
the spin at a is constrained to be J,. The corresponding 
wavefunction is 



M 



M 



M 



i<j 



U'j , (3) 



where now N is odd and M = {N — l)/2. A one-spinon 
eigenstate of Hhs is constructed by making the plane- 
wave superposition 

1 ^ 

■^„,{zi,...,zm) ^ —Y.^z*J"'^/a.{zi,...,ZM) . (4) 

Q = l 



The corresponding energy is 



1 



J.27r, 



-J—(N - —) + -( — fm{M - m) 
24^ n' 2^ N' ^ ' 



(5) 



^'m bas also a well-defined crystal momentum: q„i = 
{■k/2)N - {2n/N)(m + 1/4) (mod 27r). In terms of q^n 
the energy with respect to the ground state is E{qm) = 
{J/2)[{n/2f-ql] (modTT). 

Spinous do not lose their identity when many of them 
are present. L spinous can be thought of as a disordered 
sea with the spin at L sites constrained to be | For 
two spinous this means that the corresponding wavefunc- 
tion for a pair of localized spinous at a and [3 is given by 
(M = N/2 - 1) 

^apizi, . . . , Zm) = 



M 



M 



M 



W{Za - Z,){zp - Zj) Y[{Z, ~ Zjf Jl 



(6) 



^afi can be analytically extended to any value of Za , 
on the unit circle. As Zq,^^ are lattice sites, they are 
interpreted as locations of J,-spins. 

States with two spinous carrying well-defined crystal 
momentum are given by the lattice plane waves which 
have the expression 



V„i„(zi, . . . ,zm) - ^ 



Q,/3 



■^apizi-, . . . , Zm) 



(7) 



The total crystal momentum of ^'m„ \s q = {tt/2){N — 
2) -l- qm + Qn (mod 27r) and q™, qn are the momenta of 
each spinon. The ^'mn are an overcomplete set. A set 
of linearly independent states is constructed by taking 
only the "^mn with M > m > n > 0. Two-spinon energy 
eigenstates are linear superpositions of these: 



(8) 



where £m — n if m + n < M , £m ^ M — m otherwise. 
The coefficients a^" are |,|l^ 



(m - n + 2£) 



2^(^ + m-n+i)^ 
and the corresponding eigenvalue is 



fc-i 



(«o = 1) (9) 



En 



[E{q„,) + E{qn) 



ttJ \q„ 



qn\ 



N 



((Jm < q-n) 



(10) 



Emn is the sum of the ground-state contribution, Eqs = 
-J(7r2/24)(7V + 5/iV), and E{qm,qn), which is the two- 
spinon energy above the ground state. E{qm,qn) is the 
sum of the energies of two isolated spinous plus a negative 
interaction contribution that becomes negligibly small 
thermodynamic limit. 

The norm of can be computed by means 

of a recursive procedure, based on the operator 
ei(zi, . . . , Zm) = zi + . . . + Zm- For any wavefunction 
of the form $ x ^'gs, where $ is a symmetric polyno- 
mial, we have 



n^^GS = Egs^'^gs 



2^N^ H2 



Yz^^ 
^^dz\ 



and thus 



Zk dz. 



N -3 



(11) 
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J 2tt 



On the other hand, by standard manipulations jT^^^ , one 
can also show that 



X ^ 



GS 



iV-3 



dz 



-ei 



(12) 



^From the matrix elements of the commutator be- 
tween Tins and ei under the inner product {f\g) = 
■ ■ ■ • ■ ■ '^Af) find that, for 
the two spinon eigenstates (M = N/2 ~ 1) 



|ei|$™„) _ -{M-m+l) 



($„_i,„|$™_i,„) 2(Af-m+i) 



($™_i,„|ei|$„„) 



— (m + — n)^ 



(*m«|*mn) 2{m - u + \)m['m - Tl - \) 



($^,„_i|ei|$„ 



(13) 



IT (14) 



(^mnl^'mn) 2(n-i) 



($m.n-l|ei|$^ 



(15) 



(<I>,„,„_i|$„,„_i) 



(M - n + |)(m - n + 1)2 



2(m - n + f )(m - n + i)(M - n + 1) 



(16) 



Combining these expressions, one then finds by induction 
that 



j) r[m-n+i]r[TO-n+|] 



(*GS|*GS) 27riV(M + l)r2[m - n+ 1] 



r[TO + i]r[M 



\]T[n+\]T[M -n + 1] 



r[m+ f]r[M 



i]r[7i + i]r[Af-' 



(17) 



where ('J'gsI^'gs) = iV^^+H2M + 2)!/2*^+i |T2|. 

The definition of the wavefunction for two spinous in 
real space is now straightforward, "^ap is the state of 
two localized spinous at Za and zp. Hence, we define the 
two-spinon wavefunction, z^z^pmn{za/ zp) from 



M m 



(18) 



It is in principle possible to invert Eq.(^ and to obtain 
Pmn algebraically. However, we developed a much sim- 
pler approach, which makes use of the fact that Vta/? is 
perfectly defined for any Za, Zjs on the unit circle. Be- 
cause |$mn) is an eigenstate of Tins, one obtains 



(19) 



J,27r^2f/w d . d 9, d 

2 N y aZa OZa OZp dZj3 



2 Zq - z/3 dza azfj J 

Note the last term in this equation, which is the spinon 
interaction, is large and diverges as one power of the 
spinon separation. Upon equating Eq. ( [l9| ) to Eq. ( pO| ) we 
finally derive the differential equation 



z{l-z) 



(-m + 71+ -)z 



dz^ 



dz 



771 + 77- 

2 



Pmn — 



(21) 



The solution to Eq. (|2|) is the hypergeometric polynomial 



Pmn (-2) 



r[7n — 77+1] 
r[i]r[777 - 77 -H i] 



'^" r[fc + i]r[777 -r7-fc + i] ^ 



^-^ rffc + llr[777 — 77. - 



1] 



(22) 



In Fig.|l|we plot \pmn{za/ zp)\ VS. a — /?. The sharp max- 
imum at small spinon separation is a direct consequence 
of the strong attractive interaction between the spinous 
seen in Eq. (|2^). 

We shall now prove rigorously that this enhancement 
is responsible for the square-root singularity in the DDS. 
The susceptibility is defined by 



Xq 



M = E 



\{x\s-\^ 



GS) 



X 



{X\X) (*GS|*GS) 



2{Ex — Egs) 



{lo + irif - {Ex - EgsY ' 



(23) 



where \X) denotes an exact eigenstate of TC, Ex denotes 
its eigenvalue, and 

Sq ^Y.^O'iS^o.-^Sy) {q^2^k/N) . (24) 



However, since the act of flipping an f spin to j at site 
a is the same as creating two [ spinous on top of each 
other at site a we have by virtue of Eq. dlF 
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M m 



N 5](-l)"+>™„(l)(5(m + n - fc) 



(25) 



Thus the set of two-spinon eigenstates exhaust the ex- 
cited states coupled to '^gs by and we have 



M m 



m=On=0 



X S{m + n — k) 



2(£'mri — Egs) 



(UJ + iT^Y - {Eran - Eas f 



(26) 



This proves that the resonant enhancement is entirely 
due to the functional form of Pmn{z) shown in Fig. 1. 

The thermodynamic limit is defined as M — ^ oo, 
with m/M and n/M held constant. From general prop- 
erties of the hypergeometric functions we obtain 
p„„(l) = r[l/2] T[m - n+ l]/T[m ~n + 1/2]. Then 
approximating all the gamma functions using Stirling's 
formula and converting the sums on n and m to integrals 
over the 1-spinon Brillouin zone, we obtain the Haldane- 
Zirnbauer formula for the DSS |^] 



dqi 



dq: 



y/E{qi)E{q2 



2E{qi,q2) 



(uj + - E^{qi,q2) 



(27) 



where E{q) and E{qi, (72) are the one-spinon and the two- 
spinon energies, respectively. This may be exactly inte- 
grated over qi and q2, and the result is 



Q[^2(g) - ^] - - ^+i(g)] 

^uj - uj^i{q)^uj - uj+i{q) 



(28) 



where uj-i{q) = (J/2)q(7r - q), ij+i{q) = (J/2)(27r - 
q){q — tt) and ^2(9) — {J/2)q{2n — q). 

We see that, in the thermodynamic limit, the resonant 
enhancement in pmn turns into the square-root diver- 
gence in Xqi^) at threshold. The origin of the branch 
cut is the threshold energy for the creation of a spinon 
pair with total momentum q. The physical meaning of 
this branch cut is that the spin wave is absolutely un- 
stable versus decay into a spinon pair. Hence, no sharp 
poles, corresponding to possible low-energy spin-1 stable 
excitations, develop, but, on the contrary, the spinon-pair 
threshold is the same as the spin-wave threshold. This 



last observation points toward the main conclusion of our 
work: spinon attraction is of fundamental importance for 
understanding relevant low-energy properties of spin- 1/2 
antiferromagnets. It generates a resonant enhancement 
of the probability for two spinous to be at the same site. 
The resonant enhancement greatly increases the ampli- 
tude for a spin-1 excitation to break into a spinon pair, 
on top of an uniform two-spinon joint density of states. 
This effect is evident in the thermodynamic limit of our 
formulas, where we show that the enhancement turns into 
the branch cut in the DSS. 
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